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Abstract. Mechanisms for the generation of primordial non-Gaussian metric 
fluctuations in the context of multiple-field inflation are reviewed. As long as 
kinetic terms remain canonical, it appears that nonlinear couplings inducing 
non-gaussianitics can be split into two types. The extension of the one-field 
results to multiple degrees of freedom leads to gravity mediated couplings that 
are ubiquitous but generally modest. Multiple-field inflation offers however the 
possibility of generating non-gravity mediated coupling in isocurvature directions 
that can eventually induce large non-Gaussianities in the metric fluctuations. The 
robustness of the predictions of such models is eventually examined in view of a 
case study derived from a high-energy physics construction. 



PACS numbers: 98.80.-k, 98.65.-r, 98. 80. Bp, 98.80.Cq, 98.80.Es 

1. Introduction 

It is now clearly understood that standard single field inflation cannot produce 
significant non-Gaussianities (NG) during or immediately after the inflationary phase. 
The result obtained by Maldacena in Ref. [1 explicitly shows that standard single 
field inflation leads to no or very little primordial non-Gaussianities. And this result 
appears to be very robust, independent on the details of the model. This point is 
best illustrated by the expression of the bispectrum in the squeezed limit. Defining 
the time dependent curvature modes C(t, k) and taking advantage of the statistical 
isotropy of the universe, the power spectrum Vq of the field £ can be defined as 

(C(t,kl)C(t,k2)> = (27r) 3 <5 D irac(kl + k 2 )V ( (k 1 ,t) (1) 

and its bispectruixjj] Bq as 

(C(t, k!)C(t, k 2 )C(t, k 3 )) = (27r)^ Dirac (k 1+ k 2 +k 3 ) B c (k x , k 2 , k 3 ,i).(2) 

In the squeezed limit, i.e. when k% <C fc 2 ~ fc 3 , the bispectrum scales like 
(n s — l)V^(ki)P^(k2) where n s is the spectral index [2]. Not only are the nonlinear 
couplings naturally small - say of order unitj|§| - they are even suppressed by the 
slow-roll parameters (that ensures that n s is close to unity). 

\ In this context, both the power spectrum and the bispectrum will eventually be time independent 
at super-Hubble scales. 

§ Note however that the amount of NGs determined by a dimensionless quantity such as BfP 3 ^ 2 is 
of the order of V*(n s — 1) where V* ~ 10 — 5 is the amplitude of the metric fluctuations. 
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There are then two possible strategies to escape the limits set by Maldacena's 
results. One can modify the kinetic term by introducing higher order terms in the 
action that are not due to the potential shape. An example is provided by the Dirac- 
Born-Infeld action [3]. Such models will succeed in producing large NGs if precisely 
the kinetic term is, at the time of horizon crossing, at a non standard running point. 
That does not change however the squeezed limit case but allows large NG couplings 
for more equilateral type configurations of modes. This has been put forward as a 
powerful way for discriminating models [4] . 

Another way of evading the constraints of standard single field inflation is to 
introduce multiple scalar degrees of freedom. It can actually be argued that this 
is a natural hypothesis since it is unlikely that only one fundamental degree of 
freedom will be light (e.g. compared to the Hubble energy scale) during the epoch 
of inflation. What is more hypothetical is whether those extra degrees of freedom 
can have observational consequences. By definition, degrees of freedom that do not 
participate in the metric fluctuation, at a given time, are called isocurvature modes. 
There is no reason why the isocurvature modes should remain so all along the history 
of the universe and various mechanisms have been put forward that can lead to a 
transfer of modes, from isocurvature to adiabatic modes. 

For instance the curvaton model is based on the survival of (massive) isocurvature 
modes until late after the end of inflation that can alter the subsequent expansion 
history of the universe [5] . This is a particular case of modulated inflation [6j [7j [8] . 
Other mechanisms assume that isocurvature modes can change the end-point of 
inflation or alter the (p)-reheating effects (see [9] and contribution by A. Frolov, this 
volume). Such mechanisms can also happen in the context of hybrid inflation. It 
does not mean yet that it induces non-Gaussian metric fluctuations. That would 
happen only if there are nonlinearities in the isocurvature-curvature transfer or if 
isocurvature modes are intrinsically NG at the time of transfer. This latter situation 
is in particular advocated in Refs. [TOl [TTJ [12j [13] where isocurvature modes are shown 
to be able to develop large NG after horizon crossing. The aim of this paper is to 
show how different models that have been put forward in the literature differ in their 
mechanisms for producing NGs and how they differ in the amplitude and/or shapes 
of NGs they produce. This will be described in section 2 where it is argued that one 
can distinguished between gravity and non-gravity mediated contributions. 

The mere construction of working mechanisms cannot however be fully 
satisfactory. It is now clear that models can lead to a variety of observational 
signatures. However, whether there exist natural realizations for those models from 
high-energy physics point of view is largely open. That will be tentatively addressed 
in section 3. 



2. From single to multiple-field inflation 



The class of models we are interested in corresponds to action that takes the form 

1 



S 



dx 3 d^ 



R TOpj 



N 



(3) 



assuming that the N scalar fields $ are all minimally coupled to the metric and 
that they all have standard kinetic terms. Here, and in the following, summation 
over repeated latin indices a, 6, c, . . . from 1 to N is implicitly assumed; mp\, is the 
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reduced Planck mass, rrip 2 = 87r G. We will also assume that the spatial sections are 
Euclidean. 

As usual we will then assume that a semi classical approach can be used, e.g. 
that the fields $ a can be decomposed into an inhomogeneous part <j> a and a space 
dependent part ip a which, together with the scalar parts of the metric fluctuations, 
can be quantized (e.g. [14]). 

The zeroth order motion equation defines the field trajectory that is the time 
dependence of 4> a ■ When slow- rolling is reached it defines an adiabatic direction n in 
the field space. More precisely we can define n a as a unit vector with, 



so that n is tangential to the field trajectory during slow-roll. 

The statistical properties of the metric field and/or perturbations will then be 
obtained from an expansion of the action around the homogeneous evolution, 



where C p is of order p in the field fluctuationqJj|. 
2.1. The equations 

From a practical point of view the choice of gauge can be crucial. Here we chose a 
spatially flat slicing gauge and derive the motion equation using the ADM formalism. 
This technique is very efficient and has been presented (and used) in Ref. [1] in this 
context. We refer the reader to this paper for a detailed presentation of the method. 
In the ADM formalism the metric is written, 




(4) 




(5) 



ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx :i + N j dt). 

and we make the gauge choice so that 

hij = a 2 (t)5ij, 

introducing a spatially flat slicing gauge. 

The resulting second order action is then given by, 

£ 2 = - Vab <PatPb + fil - -o {difaf 

n 



(6) 



(7) 




(8) 



The motion equation of the field ip a can then be read out of the action, 



(f a + 'dHipa ?A(fi a = S ab (f b 

ft 



(9) 



a 



with 




(10) 



(11) 



|| Up to a total derivative term, C\ vanishes about the homogeneous trajectory. 
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In the slow-roll regime, it can then be written, 

C = —7= n a tp a , (12) 

TOpi. \J it 

where e is the (first) slow-roll parameter, 2e = (m P1 n a V a /V) 2 . This set of equations 
provides the necessary ingredient to compute the shape of the power spectrum however 
complicated the potential might be. 

To obtain the induced amplitude of the NGs it is however necessary to extend 
the calculations to higher order terms. Both the relations © and (jlll) should then be 
extended to higher order terms. These calculations were initially carried by Maldacena 
in Ref. [JJ for a single field. 

The relation between the curvature perturbation and the field fluctuations is in 
general intricate. The SN formalism offers however an insightful way of representing 
those couplings ([15l [16] and see Ref. [17] in the context of nonlinear expansions). 
Indeed scales of interest being super-Hubble, the curvature can be identified as 
fluctuation of the number of efoldings from one field trajectory to another. This 
number can in turn be expanded as a function of the field fluctuations, 

TT^afc <Pa(t*)<Pb(t*) + ■ ■ ■ (13) 



SN(t) = N, a ip a (U 



where 



N a = 



dN(t) 

dip a (t*) 



and N. ab = 



d 2 N(t) 



(14) 



d<p a (U)d(pb{t*) ' 

Here is a time for which all modes are evolving at super-Hubble scales. This is 
usually taken as the time of horizon crossing (as in Refs. [TBI HI HO]) but this is not 
necessarily so and it can be taken latefl The presence of a second order term ensures 
that the resulting metric fluctuations are not Gaussian. This is however not the only 
possible contribution. Other possibilities come from the fact that the fields themselves 
may not be Gaussian distributed to start with at time t*. The unambiguous way to 
derive the statistical properties of those fields is to start with the third order action 
and, using for instance for the "in-in" formalism from Schwinger and Keldysh (see 
Refs. [UJ [25] and [23] for a presentation and use in a cosmological context), derive 
the properties of the field. 

In case of multiple-field inflation the result obtained by Maldacena for single field 
can be readily extended. The third order action is given by (using the same gauge as 
before) , 



1 
6 



4Hm%, h 



8Hn4 h 



Vbc tfb<p c + fib<Pb + -^ditfbdttpb 
a z 



2H 2 



{<t>aVa)' 2 

AH 



— 4>bVb + d 2 \ 



AHm% x 



4aVa [didjxdidjX - d 2 xd 2 x] , 



(15) 



% That amounts to compute the nonlinear super-Hubble evolution of the isocurvature fields either 
separately, from the action expression, or to include it in the JV a t coefficient as illustrated below. 
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where 

d\ = 2H^T^- (16) 

In Eqs. ([T5]) and ([T6"]) the different terms are ordered in terms of the power of 
the Planck mass. They can alternatively be ordered with respect to the slow-roll 
parameters. For single field inflation such a decomposition is unambiguous as described 
in [TJ . For multiple-field models the decomposition depends on whether the derivatives 
of the potential have the same order of magnitudes in all field directions. If it is so, 
the third order action is dominated by the last three term of the second line: they are 
all of order \/eH /mpi.. That corresponds to the expression of the action as derived in 
[24] . On the other hand, if only the adiabatic direction has small derivatives then the 
expansion should be made with more care. In particular, the first term of the second 
line can be of the same order of the other terms of that line; and the very first term can 
dominate all others inducing non-gravity mediated couplings. Indeed, whereas field 
fluctuations in the adiabatic direction are intimately coupled to metric fluctuation, 
and are therefore constrained by slow-roll conditions, this is not the case in transverse 
directions (provided its effective mass remains below H, see discussion in Ref. [12]). 
and field can develop arbitrarily large couplings. To give another insights into the 
separation of those terms, one can note that in the de Sitter limit - which implies that 
all time derivatives of background quantities vanish - the only non-vanishing term 
of £3 is the first one. It simply describes the self-coupling of fields in an arbitrary 
expanding background. 



2.2. A geometrical description 

We have now all the necessary ingredients to present the various mechanisms at play 
in this context as is illustrated in Fig. Q] Those sketches describe field trajectories 
assuming one isocurvature degree of freedom. The right part of the panels correspond 
to sub-Hubble dynamics where quantum fluctuations are at play. Non-Gausssianities 
can be induced at that stage as described by the "in-in" formalism. The formula (fT"7f 
below gives the expression of the field bispectrum but with the assumption mentioned 
before, that is assuming slow-roll conditions are met in all field directions. It docs not 
therefore encompass all possibilities. 

The transfer from isocurvature to adiabatic modes can take place as soon as 
modes cross the horizon or much later on. In this first case, a full treatment of 
the metric perturbation with the help of the "in-in" formalism is required. In the 
latter case further nonlinearities are built during the super-Hubble evolution of the 
fields. Those evolutions can be taken into account through the SN formalism (which 
corresponds to a classical evolution of the fields) or in the "in-in" formalism applied 
to later time. Both descriptions are equivalent and should (see Refs [25] and [26] for 
a general discussion on that) give the same result. 

It is to be noted then that if the super-Hubble evolution of the field is the 
dominating process, then the bispectrum is going to be local. The later case can 
be further extended to situations where mode transfer takes place at the very end 
of inflation. This is the case in particular for extensions of the hybrid models (see 
Refs. [Ml E [271 ES [TO]). Tw0 cases though should be distinguished: when the 
transfer is itself nonlinear (that would correspond to critical region with a strongly 
curved boundary on the bottom panel) or the transfer is linear (the boundary is 
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Figure 1. Sketch of the sub- to super-Hubble evolution of the field fluctuations 
in isocurvature and adiabatic directions and the mechanisms that can lead to 
significant non-Gaussian metric perturbations. In the top panel, the transfer of 
modes takes place during the inflationary period while in the bottom at the end 
of inflation assuming it corresponds to the critical region in grey. 



essentially a straight line) in which case the non-Gaussianities are those imprinted in 
the isocurvature modes. 

This latter picture is the one put forward in p~3j [27] and advocated as a very 
efficient mechanism for producing non-Gaussianities since it can take advantage of 
non-gravity mediated mode couplings. 



2.3. Gravity mediated couplings 

When only the leading order in slow-roll parameters is included, and assuming the 
potential derivatives are of the same order in all directions, one can derive the field 
three-point correlation functions (e.g. [24 | 129 ] ITS]). 

(^a(kiV 6 (k 2 V c (k 3 )) = (27r) 3 <5 D irac(kt) B a6c (k 1 ,k 2 ,k 3 ) 



= (2vr) 3 fc irac (kt) 



A? 



perms 



M(k 1 ,k 2 ,h) 



4>aS, 



be • 



with 



M(ki,k 2 ,k 3 ) = -k\kl 



h 2 h 2 



h 2 h 2 



+ + — ■ 



(17) 



(18) 
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where k t = k! + k 2 + k 3 , and it is assumed that all the k t are of the same order of 
magnitude so that they cross the Hubble radius approximately at the same time. As 
a result A 2 = k 3 P ip (k) is, at leading order in slow-roll parameter, the same for all 
three k^ at horizon crossing. The sum that appears in (|17[) is over all simultaneous 
rearrangements of the indices a, b, and c, and of the momenta kx, and £3 in A4, 
such that the relative position of the hi is respected. The form presented here for the 
function A4 is the one given in Ref. |18j rewritten from the expression given in Ref. 

It can be noticed that the resulting field bispectra scale like 1/^pj and are 
proportional to the field time derivative. As a result, when re-expressed in terms of 
the slow-roll parameter, they are typically of the order of 

B«6c(ki,k 2 ,ka)~-^-P£ (19) 
mpi. 

which implies that the bispectrum of the metric fluctuations is of the order of 

B c ~eV 2 (20) 

not taking into account the subsequent growth of mode couplings. What we recover 
here is qualitatively nothing but the result obtained by Maldacena for single field 
inflation in Ref. pQ. 

2.4- Bispectra from non-gravity mediated couplings 

By construction the previous section ignored the non-gravity mediated terms. Here 
we rather focus on the other limit case, that is we assume that fields in the transverse 
direction develop large NG through intrinsic couplings (e.g. first term of Eq. (fT5l) ). 
So let us consider N — 1 isocurvature modes Xi (which is therefore a N — 1-dimension 
subspace of the initial iV-dimension field space) and we keep here the only cubic term 
of the potential (mainly to be on the same footing with the previous results), 

d 3 V(iJ,Xi, -,Xn-x) = Xi JK (v) ^ 
dxiXJXK 3! 

It is then possible to show explicitly that the x fields develop a non-zero bispectra at 
horizon crossing and after. The amplitude of the coupling is determined by the values 
of Xijk that a priori have no relation with the slow-roll parameters. It can be shown 
that it takes the form (see Refs. [12] and j3U] for the general case), 

Bzjx-(ki,k 2) k 3 ) - -vi j K [P{k 1 )V(k 2 ) + sym.] , (22) 

with, 

vijk = J dr/aV)A/^(?/) J _ ^T y (23) 

which, in case of a de Sitter background and when the coefficients Xijk are constant 
in time, reduces to 

(n c - jy. 

3H 2 

This expression is valid for sufficiently late time after horizon crossing. In case of the 
de Sitter case that simply means that the number of cfoldings N c — N* since horizon 
crossing for the scales of interest is large (i.e. N c — ^> 1). The results that are 
presented here correspond to cases where the super-Hubble evolution domimates over 



vijk = Xijk " ° olT2 * ■ ( 24 ) 
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coupling effects at horizon crossing time. This can be shown explicitly using the "in- 
in" formalism (an early derivation of the three-point correlation function of a test field 
in a de Sitter background can be found in [3T], the four-point correlation in [52]). 

An important consequence of this approximation is that the bispectra take a local 
form, e.g. the dependence of the bispectrum is simply that of a product of spectra. 

What is then the amplitude of the couplings? In case of a cubic potential one 
needs the effective mass of the isocurvature mode to be small compared to the Hubble 
constant. That imposes that the elements of the matrix 1/ ' H^jXijk are small 
compared to unity, the (pj being all of the order of -H0. 

That implies that the bispectra can be up to the order of 

£ x (k 1; k 2 ,k 3 ) < (^-^*) [p^VM + sym.] (25) 

This is to be compared to Eq. p^|) . Here, the resulting amplitude is not suppressed 
by the slow- roll parameter anymore; it is also a factor mp\ ,/H larger! That shows 
how strong non-gravity mediated couplings can be. There is however a missing 
ingredient: this mechanism can only be efficient if a transfer of modes occurs during 
the inflationary period of the universe, or at the very least at the end. 

3. Multibrid inflation, from toy models to SUSY models 

What the previous analysis has shown is that multiple-field inflation offers plethoras of 
mechanisms to play with. In the absence of guideline from high-energy physics it seems 
therefore difficult to extract robust common features. Attempts to overcome this issue 
with the help of high-energy motivated constructions have all, in this context, lead to 
inflationary models that are hybrid and or extensions of hybrid type inflation. Hybrid 
inflation (see Ref. [3 5) ) is by itself an appealing model; its extensions in the context 
of multiple-field inflation can lead to interesting phenomenological properties we wish 
to describe in this section. Let us start with a simple model we will later motivate 
from super-symmetric (susy) theories. 

3.1. A toy model 

Let us consider a simple extension of hybrid inflation with one field (this model has 
been proposed in Ref. [IS]), 4>, the inflaton; a light scalar, x> that will give birth 
to isocurvature fluctuations and a third field, a, which is coupled to the two others 
so that the end of inflation is triggered when a undergoes a phase transition. To 
make the model robust in particular with fine tuning in parameters and in the initial 
conditions, we introduce a quartic coupling in the isocurvature direction so that the 
potential reads 

V = V{4>) + ^x 4 + ~ (a 2 - a 2 a ) 2 + fa 2 (0cos0 + X sin0) 2 . (26) 

Here Co is the vacuum expectation value (vev) reached by a after the phase transition, 
9 is the mixing angle between <\> and \ in their coupling to a, v 2 is the amplitude of the 

+ Another way of arriving to a similar amplitude is to start with a (positive) quartic potential of the 
shape Ax 4 - In this case A (which is then dimensionless) ought simply to be less than unity. Finite 
volume effect though gives a non-zero value of \ corresponding to the one it acquires at the size of 
our observable universe. Using the Fokker-Planck approach described in 33 it can be shown that it 
takes typically a value of about H/X 1 ' 4 1341 . It leads to a non-zero cubic term of similar amplitude. 
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Figure 2. Expected one-point probability distribution distribution of the 
isocurvature perturbation in case of a quartic potential (solid lines). The plot 
corresponds to the case v 2 (N c — 7V*)/9 = 1. Finite volume effects can shift the 
Hubble size value of \ to a non-zero value XHorizon for our observable universe 
leading to a skewed distribution of \- The three solid lines correspond to 
XHorizon/^ = 1, 0.5 and 0. The dashed line is a Gaussian distribution of the 
same width. See Rcf. 34 for details. 



self-coupling. We assume waterfall conditions so that inflation effectively ends when 
the effective mass of a vanishes, that is when 



The value of < 



g(4> cos 9 + x sin 9) — 2^o-q = 0. 
at the end of inflation is thus 



fend 



cos( 



(27) 



(28) 



For 4> > </>end, cr — and the two fields </> and \ evolve independently: cf> drives the 
inflation while x develops non-Gaussianity. The amount of non-Gaussianity of x then 
depends only on v and on the total number of e-foldings between horizon crossing 
and the end of inflation. 

When 9 is non-zero, fluctuations of x induce metric fluctuations because they 
change the time at which the phase transition occurs. Thus the x~ m duced metric 
fluctuations read (assuming H is basically constant during the inflationary period), 

, TTC 3H 2 sin 9 

(~H6t cnd ~-— -x- 29 

V<p cos 9 

Here the transfer of modes is linear and the induced metric fluctuations inherit their 
non-Gaussian properties from those of the isocurvature modes. 

It is to be noted that this model does not require any fine tuning to work, neither 
in the parameter space (the only requirement is that v is less than unity), nor in the 
initial conditions (if initial conditions are such that x lS large, this field is bound to 
roll down near its minimum at a vev of about H). It is actually to be noted that, 
despite the fact that the isocurvature potential is symmetric, finite volume effects 
induce generically a non zero bispectrum but the amplitude of which depends on 
the peculiar (patch of the) universe we live in. This mechanism has been described 
more quantitatively in 34 . The picture we eventually arrive at is that the metric 
fluctuation is the superposition of a Gaussian field and a non-Gaussian field the latter 



Gravity and non-gravity mediated couplings in multiple- field inflation 



10 



being obtained essentially from a local nonlinear transformation^! of another Gaussian 
field of the same variance, 



with 



Glxo] 



3H 2 



Xo 



3H 2 



V, 



tan(9£[x ] 



g H 2 



-1/2 



(30) 



(31) 



where xo is a Gaussian field of variance H* (but not necessarily of zero mean). Fig. 
[2] illustrates the expected shapes of the one-point probability distribution function 
(PDF) of x that one obtains from such a transform. 



3.2. Protecting masses in multiple- field inflation 

From a high-energy physics point of view, the introduction of such multiple degrees 
of freedom raises a number of questions. Which constructions are more natural? Are 
such models robust in a high-energy physics context? In particular the protection of 
the masses and coupling constants of the scalar degrees of freedom we introduced is 
critical for the viability of the models. If large-scale divergences can be taken care of, 
mainly from horizon effects or invoking the semi-classical approach described by the 
Fokkcr-Planck equation, [33] , masses and coupling constant are a priori unprotected 
from small-scale divergences, a problem encountered in expanding backgrounds as 
well as in Minkowski space-time. In case of a quartic potential, it is always possible 
to invoke renormalization properties but it then requires fine-tuning. To avoid it it 
has been shown in Ref. |36j that if scalar degrees of freedom were imbedded in a 
super-symmetric multiplet, their power spectrum - and as a consequence their masses 
- are protected from radiative corrections. Construction in a susy context is then a 
natural playground for the derivation of viable models and one reason that makes 
models such as D-term or i^-term inflation described respectively in Ref. [37] and in 
Ref. [38] attractive. 



3.3. A viable model in the context of D -term super- symmetry 

It goes beyond the scope of this short review to describe the mechanisms at play in 
the construction of the D— term inflation (see for instance Ref. [55] )• It is based of 
the up-lifting of flat direction due to (weak) radiative corrections when the current 
vacuum breaks a supersymmetric Lagrangian. In case of D-ierm. infation, the whole 
structure derives from the expression of the superpotential that in this case takes the 
form 

W = \S~4><j) (32) 

where S, <f> and <j> are neutral or charged fields and A a dimensionless coupling constant. 
The resulting potential is 

V = V^ loop (S) + A 2 |S| 2 \W + ^ (-M 2 + C) 2 (33) 

* Note that strictly speaking it can be shown that the tree order correlations of the induced \ field 
are identical by those induce by such a transform (see Refs. 1121 1301 ). If radiative corrections were 
to be taken into account the above statement would not be correct. 
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Figure 3. Exclusion diagrams for parameters v and % for 8 = it/A (left panel) 
and for 9 = 0.1 (right panel). The gray areas are the exclusion zones, dark gray 
corresponding to the constraints derived from WMAP and medium and light gray 
for the constraints expected from Planck. The straight line boundaries correspond 
to bispcctrum observations (of slope —2), the other case to the trispectrum. The 
long dashed is the location where the terms contributing to the trispectrum cancel. 
We adopted the results of 141 1 on the upperbounds the Planck mission is expected 
to provide, /nl = 5 and t>jl = 560. The dotted lines are the locations where \ 
is equal to its expected 1-cr fluctuation. 



where g is a gauge field coupling constant, £ 1 ' 2 if the energy scale provided by the 
Fayet-Iliopoulos term and which depends on the (complex scalar) fields S and <\>. Here 
Vi-ioop is the contribution due to the radiative corrections (computed from the formula 
of Coleman and Weinberg, 40 J that ensures a slow rolling of S towards its minimum. 
One recognizes here a hybrid model. 

But then, following [27] , nothing prevents the introduction of multiple light fields 
Si coupled to the same charged U(l) fields, 

w = E t s ? + A (E a ^J ^ ( 34 ) 

where Vi and ai are dimensionless parameters. Obviously if a peculiar self-coupling 
parameter v\ is small enough, the corresponding Si field can participate in the inflaton 
(depending on its initial vev). The corresponding upper bound for v\ for such a 
possibility to occur comes from the fact that when the vev of Si is below the Planck 
scale, the contribution of the quartic potential it induces should be negligible against 
the radiative correction terms. It leads to the constraint, 

vl « A 4 . (35) 

The fields for which Vi is above this bound will rapidly roll towards the origin but 
they still can develop significant supcr-Hubble fluctuations as long as Vi is smaller 
than unity. 

The resulting expression of the potential corresponds to the following effective 
form, 

V = Vl-ioop + ^ 2 2 |5 2 | 4 + A 2 |cos6»5i + sin<9S 2 | 2 1^| 2 

+ ^(-H 2 +e) 2 06) 



involving the (complex scalar) fields Si, 52 and <f>. It corresponds to the model ([26| 
with <p = Si, x = ^2 with the difference that the fields are complex (and thus carry 
two degrees of freedom each) . It changes the resulting phenomenology only in details 
(for instance on the expected amplitude of the trispectrum, see Ref. [27]). 
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The entire non-Gaussian properties of such a model are determined by three 
parameters that are, 

• the coupling constant v that sets the amplitude of the non-Gaussianities that 
can be reached in the isocurvature mode. For perturbation theory to be valid, 
v 2 should be less than unity. Actually we found that the combination N e v 2 , 
where N e is the number of efoldings since horizon crossing, ought to be small. 
This condition ensures that the field \ is light compared to the Hubble constant, 
its effective mass is v 2 H 2 , and that its evolution is perturbative all through its 
super-Hubble evolution. 

• the mixing angle 9 that determines the fraction of isocurvature modes that 
eventually enters the metric perturbations. This angle can be anywhere between 
and tt/2. 

• the super-Hubble value of x, x, is an priori undetermined value which depends on 
the patch of the universe we live in. Within some hypothesis, its cosmic PDF can 
be determined and therefore its expected values. Such derivation is done with 
the help of the Fokker-Planck equation (see Ref. [27J in this multi-dimension 
context). The value \ determines in particular the amplitude of the expected 
bispectrum (it is zero if x is set to 0) and the relative contribution of the terms 
contributing to the trispectrum. 

On Fig. [31 one can see the current and expected constraints on parameter space 
that can be put on such a model in the v — \ plane for two different values of 6. The 
upper right parts of the diagram are a priori excluded. It can be seen that bispectrum 
and trispectrum play a complementary role in constraining these models. 

4. Conclusions 

The exploration of the various types of coupling terms that appear in the action leads 
to distinguish gravity from non-gravity mediated couplings. In one-field inflation, 
because the field fluctuations and the metric fluctuation are locked together, only 
the former can be found. In multiple-field inflation however this is not necessarily so 
and it opens the possibility of having a richer phenomenology. Whereas the gravity 
mediated couplings are ubiquitous but induce only modest effects, the non-gravity 
mediated couplings can be very efficient, although nothing ensures that they are 
generically at play. Giving the freedom one has in building potentials it is however 
certainly possible to design models exhibiting any kind of scale and geometrical 
dependences in the bispectrum (this has been tentatively explored in Ref. 42]), but 
actual constructions motivated by high-energy physics point to hybrid inflation type 
models where couplings are eventually of local type. The extended D-term susy model 
presented here is an example of such a construction. A few lessons can be drawn from 
its analysis, 

• in such a susy framework masses and coupling constants are naturally protected. 
This is at the expense of the doubling of the number of scalar degrees of freedom; 

• this extension of D-term inflation leads to hybrid models where the end line, 
the critical line where the inflationary period terminates, is linear in all field 
directions. This is at variance with the construction proposed in [2H] where NGs 
originate from a nonlinear critical line; 
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• in such a model there is no need for specific fine tunings, neither in the initial 
conditions nor in the parameter space. It is also to be noted that such models 
induce a dumping of the rare event tails - on both sides. 

• the transfer of modes, from isocurvature to adiabatic direction, takes place at the 
time the inflation stops. It leads to local types for bispectrum, and trispectrum, 
shapes. 

• calculations were carried here at leading order in perturbation theory. The 
conditions for such calculations to be valid during the sub-Hubble evolution is that 
the coupling constant v 2 is less than unity. Once the evolution is super-Hubblc 
however the amplitude of the fluctuation couplings is driven by v 2 {N ~ — N*) which 
ought then to be small. There is therefore a regime where perturbation theory 
can be used during sub- Hubble evolution but not during the whole super-Hubblc 
evolution. That could lead to a new phenomenology; e.g. to nonlocal effects in 
the shape of the high order spectra. 

The model at hand is therefore rather sound with not so much freedom in its 
predictions. But this is by no means exhaustive. Other scenarios are certainly possible. 
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